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Abstract. There are various ways of measuring the slope of the upper end of the 
IMF. Arguably the most direct of these is to place stars on the H-R diagram and compare 
their positions with stellar evolutionary models. Even so, the masses one infers from 
this depend upon the exact methodology used. I briefly discuss some of the caveats and 
go through a brief error analysis. I conclude that the current data suggest that the IMF 
slopes are the same to within the errors. Similarly the determination of the upper mass 
"limit" is dependent upon how well one ca n determine the masse s of the most massive 
stars within a cluster. The recent finding bv lCrowther et al.l (1201 Ot) invalidates the claim 
that there is a 15OM upper limit to the IMF, but this is really not surprising given the 
weakness of the previous evidence. 



1. Introduction 

The goal of this talk is to summarize what I think we know (and don't know!) about 
the upper end of the IMF and the upper mass limit, with an emphasis on how we know 
what we know. The latter allows us to understand how well we know what we think 
we know. When the organizer's invited me to this conference I got hooked by their 
statement that we would "critically re-evaluate the ensemble of accumulated evidence 
constraining the upper end of the IMF." I come here as observer wishing to remind my 
colleagues of the uncertainties and caveats that are sometimes overlooked. Or to put it 
in NPR's Car Talk parlance, I'm here to administer an observer's dope slap, where a 
dope slap is "an attention getting device." 

2. IMF Slopes of Clusters Containing Massive Stars 

Let me begin with the most "direct" way of measuring the IMF: counting stars in the 

Hertzsprung-Russell diagram (HRD). 

Table Q] lists some IMF slope determinations fr om iMassey ( 1998), updated with 



more recent values for R136 (Massey & Hunt er! 19981) and h and^f Persei ( Slesnick et al. 



2002). Included are the ages inferred for the clusters and the masses of the highest mass 
present. These numbers were all determined in a uniform way, and there are certainly 
differences between the IMF slopes T's that are larger than the quoted errors. But seri- 
ously, where do these numbers come from, and how much can one believe them? Do 
these numbers demonstrate that the IMF slope indeed varies at the upper end? 

The errors quoted on the IMF slopes T in the above table are strictly the fitting 
uncertainties to the number of stars as a function of mass. They do not take into account 
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Table 1 . Ages, Upper Mass, and IMF Slopes 



Region 


Age 
|MyrJ 


M up 
r A/f ~\ 


1 
1 


Ker. 








SMC 




NGC 346 


2-4 


70 


-1.3+0.1 


Massev et al. (1989a) 








LMC 


LH6 




85 


-1.7 + 0.4 


Oev ( 19961 


LH 9 


1-5 


55 


-1.4 + 0.2 


Parker et al. (1992) 


LH 10 


0-3 


90 


-1.1 ±0.1 


Parker etal. (1992) 


LH38 




85: 


-1.6 + 0.2 


Oev (1996) 


LH 47/48 


2-3 


50 


-1.3+0.2 


Oev & Massev (1995) 


LH58 


2-4 


50 


-1.4 + 0.2 


Garmanv etal. (1994) 


LH73 




65: 


-1.3+0.4 


Oev (1996) 


LH83 




50: 


-1.3+0.5 


Oev (1996) 


LH 114 




50: 


-1.0 + 0.1 


Oev (1996) 


LH 117/118 


1-3 


100 


-1.6 + 0.2 


Massev etal. (1989b) 


R136 


1-2 


>120 


-1.3+0.2 


Massev & Hunter (1998) 








Milky Way 


NGC 6823 


2-7 


40 


-1.3+0.4 


Massev et al. (1995) 


NGC 6871 


2-5 


40 


-0.9 + 0.4 


Massev et al. (1995) 


NGC 6913 


4-6 


40 


-1.1 ± 0.6 


Massev et al. (1995) 


Berkeley 86 


2-3 


40 


-1.7 + 0.4 


Massev et al. (1995) 


NGC 7235 


(6-11) 


15 


(-2.0) 


Massev et al. (1995) 


NGC 7380 


2 


65 


-1.7 + 0.3 


Massev etal. (1995) 


Cep OB5 


2-4 


30 


-2.1 +0.6 


Massev et al. (1995) 


IC 1805 


1-3 


100 


-1.3+0.2 


Massev et al. (1995) 


NGC 1893 


2-3 


65 


-1.6 + 0.3 


Massev et al. (1995) 


NGC 2244 


1-3 


70 


-0.8 + 0.3 


Massev et al. (1995) 


NGC 6611 


1-5 


75 


-0.7 ± 0.2 


Hillenbrand etal. (1993) 


CygOB2 


1-4 


110 


-0.9 ± 0.2 


Massev & Thompson (1991) 


Tr 14/16 


0-3 


>120 


-1.0 + 0.2 


Massev & Johnson (1993) 


h and x Per 


12.8 


35 


-1.3+0.2 


Slesnick et al. (2002) 



any of the other errors, such as the determination of the masses, or for that matter, the 
degree of coevality that is needed to equate the slope of the present day mass function 
with that of the initial mass function. Let us consider at some length just one of these 
issues, namely the determination of the masses. 

2.1. Determining the Masses 

To count the number of stars as a function of mass, one needs to somehow get the 
masses. To do this, we must find the star's bolometric luminosity and then rely upon 
stellar evolutionary models to convert the luminosity to mass. Generally one begins 
by obtaining UBV photometry of the stars. These are used to identify the hot, lumi- 
nous stars, usually by the Johnson reddening-free Q index (U - B) - 0.72(fi - V). One 
then takes spectra of as many stars as practical, and uses each star's spectral type to 
convert to effective temperature T e ff . The effective temperatures yield the bolometric 
correction as a function of spectral subtype. Comparison of a star's observed colors 
with that expected based on the spectral type and luminosity class also yields the color 
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excess and hence the correction of V for reddening. In the cases where the distances 
are poorly known (as in Galactic clusters) the spectral types can also be used to de- 
termine a "spectroscopic parallax" to the region, i.e., using the observed de-reddened 
magnitude and the absolute visual magnitude expected from each star's spectral type. 
For clusters older than about 10 Myr it is more accurate to use main-sequence fitting 
from the photometry itself t o determine the distance modulus (as we did for h and x 
Per; see lSlesnick et al.ll2002l) . 

So, are spectra really necessary? Couldn't we just go from O to T e ff and be done 
with it? Let's consider the associated errors, following iMasseyl ( 1998 ). First, assume 
that the mass-luminosity relationship goes something like L ~ m 3 , in accord with the 
solar-metalllity Geneva evolution models for stars with masses > 60M Q . (Below this 
mass the exponent is larger.) Then Alogm ~ 0.3AlogL. If we express the luminosity 
in terms of the bolometric magnitude M\, \, then A log m ~ O.lAMboi- 

The uncertainty in M^ \ will be a combination of the error in the reddening, the er- 
ror in the distance, and the error in the bolometric correction BC. The latter is related to 
the effective temperature 7V ff as B C = -6.90 log T eS + 27.99, so ABC = -6.9Alog T eS , 
according to lMassey et all (120051) . 

Let us now see how the uncertainties in log m stack up. The typical uncertainty in 
the distance is 0. 1 mag, and so that is one source of uncertainty on log Mb i, but of all 
of the errors that one will be systematic for stars of all masses in a cluster and will have 
the least overall effect on the derived IMF slope, but it will affect the derived masses 
so let's include it. The reddening might be uncertain by a similar amount, and in this 
case it can be either systematic or vary star-by-star. If we are using spectral types to get 
log r e ff, then an uncertainty of 1 spectral type corresponds to an uncertainty of about 
0.01 in log T e g and hence 0.1 in log Mbol- Thus all three terms contribute similarly, and 
the total uncertainty in M^ \ is about 0.2 mag. The uncertainty in the mass (assuming 
the evolutionary tracks are gospel) is then 0.02 dex, which is really very good! (We will 
talk about systematics below.) 

If instead we had depended upon the photometry to derive a mass, we would have 
to have related a reddening-free index (such as Q) to log T e s using stellar atmospheres. 
For hot, massive stars this is a very frustrating endeavor, as the effective temperatures of 
these stars are so high that there is very little difference in the c olors of a 50,00 K star 
and a 40,000 K star, even though the spectra are quite different. IMasseyl ( 1998b derives 
values for Q and other reddening-free indices, including ones derived from spacecraft 
UV-filters. I found that the standard UBV Q would vary by only 0.007 mag over the 
same temperature range. In general, for hot stars, A log T e g ~ 5AQ. So, a (reasonable) 
error of 0.05 mag in Q would translate to an uncertainty of 0.25 dex in log r e ff, and (via 
the equation for the BC), about 1.7 mag in log Mt, i, or 0.2 dex in log m, about lOx the 
error from using spectral types. You can't really tell a 100M Q from a 60M Q star using 
photometry alone. (I'll briefly note that the UV filters fare a bit better but not much.) 
This is all better illustrated in Figure [T] where the relative errors of placing hot, massive 
stars in the H-R diagram is shown. 

How would these errors translate into errors in the IMF slope T? If the errors were 
somehow systematic in log m they would actually have no effect. But far more likely 
there would be biases introduced as a function of mass, as the hottest stars (where the 
dependence upon photometry is even less) are also the most massive. 

We can demonstrate this by just comparing how different researchers using dif- 
ferent methods agree on the IMF slope. Consider the value for LH 58, a beautiful OB 
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Figure 1 . On the left is shown the errors associated with placing stars on the H-R 
diagram using their spectral types to determine their masses assuming an uncertainty 
of one spectral subtype. On the right is shown the errors associated with placing the 
same stars based upon their photometry, assuming an uncertainty of 0.05 in Q (i.e., 
roughly 0.02 mag error each in (U-B) and (B-V). 



assoc iation in the LMC. It was a nalyzed by ourse lves (IGarmany et all 1 19941 : iMassey 
1998) using s pectroscopy, and by Hill et al.1 dl998l) using their own photometry bu t no 
spectroscopy. IGarmany et al. (1994) derive r = -1.7 ± 0.3, while iMassevI (Il998l) de- 
rives T = -1.4 ± 0.2 using the identi cal data, but tre ating the reddening jus t a little 
bit dif ferently. Using only photometry, iHill et all (1 19981) find T = -2.5 ± 0.3. iMassev 



(1 19981) notes that reanalyzing the cluster adopting the Hill et al. ( 1998) conversion from 
photometry to effective temperatures, but using the IGarmany et all (119941) photometry, 
leads to T = -2.0. The conclusion I take away from this is that reasonable people using 
reasonable data may still derive IMF slopes that disagree by a significant amount! If 
one is going to look for variability of the upper end of the IMF this way one has to 
analyze the data consistently. Figure [2] shows the IMF slopes from Table 1, where the 
data have been treated as much the same as possible. It's hard to come away from this 
thinking there's evidence for a variable IMF slope. 

One should recall that while these values for T have been found consistently, that 
doesn't mean that they are right in some absolute sense. For one thing, binarity has 
been ignored. The relative values are still right if the binary frequency and mass ratio 
distributions are the same from region to region. But the actual values shouldn't be. 
Furthermore, our conversions have relied upon a set of conversions from spectral types 
to effective temperatures. That s cale is probably uncertain by 10% in an absolute sense 
(IContilll998l : lMassev et al.ll2005h . 

I should also note that the way we went about determ ining the IMF slope once we 
had the masses was a bit crude. As [Bastian et all d2010l) notes, one could do perhaps 
better by a more sophisticated analysis, comparing a chi-square goodness of fit with 
models of the cumulative distributions drawn from various IMF slopes. 

What about determining the IMF slope from just a luminosity function? Well, that 
turns out to be incredibly insensitive to the value of T, at least in the case of mixed-age 
populations. Consider the case of looking at massive stars in a nearby galaxy. One 
can't use all of the stars in a color-magnitude diagram as the vast majority of these will 
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Figure 2. IMF slopes determined in as consistent a manner as possible are shown 
for OB associations and clusters in the SMC, LMC, and Milky Way. The errors 
are undoubtably larger for the Milky Way regions due to th e addit i onal reddening 
and un certa inties in distances. The d ata come from Massey ( 1998), Slesnick et al. 
(2002), and lMassev & Hunted (fl998l) . and cover a range of 4 in metallicity and fac- 
tors of several hundred in stellar density. A Salpeter IMF slope of T = -1.35 is 
shown for comparison. 



be do minated by yellow foreground stars (see, for example. iMassewOlOHDrout et al. 
120091) . One can restrict the sample to the bluest stars only, say, B - V < 0. But the 
sensitivity to the IMF slope is almost non-existent, as shown in Figure [3] This works a 
lot better for stars that are coeval, but without first constructing an H-R diagram it's not 
clear how one knows a priori that a region is coeval. 



3. The Upper Mass "Limit" 



Massevl (l998) wrote that the highest mass stars identified in OB clusters and associ- 



ations were consistent with what one expected by extrapolating the IMF slopes of the 
lower mass stars : that generally any " truncation" in the mass function was simply due to 



age. In contrast, lOey & Clarke! d2005r) argue that by combining data from many regions 
(mostly those in Table 1) that one should expect purely by chance to have seen higher 
mass stars than what have been observed so far, and argue that stars more massive than 
150M o must not be able to form. 

Let us think about what the requirements would have to be to find a (say) 300M o . 
First, my own expectation would be that we would find it in a rich association, not a 
sparse one. Not everyone would agree with that point, as there is another school of 



6 



P. Massey 




Figure 3. The luminosity functions are amazin gly insensitive t o the IMF slopes, 
as demonstrated for a mixed-age population. From Massey] (119981) . 



thought that somehow the probability of forming a star of mass m is independent of 
physical conditio ns. I confess I've n ever understood that argument, but that doesn't 
mean it's wrong. iBastian et all (2010) presents nice arguments relating the mass of the 
highest mass star seen to cluster size. But, regardless of whether we all agree on this 
point or not, we probably can all agree that if were hoping to find a very high mass star 
in a cluster, we would probably do better to look at clusters that were sufficiently young 
to still have such a star among the living. 

How young is young enough? Here we are on somew hat shaky grounds, as the 
best yetted models are those of the Geneva group (for instance. iMevnet & Maede i]|2005l 
2003| : ICharbonnel et al.ll988l : ISchaerer et alJl993l : ISchaller et all 1 9921 : iMaeder & Mevnetl 
19891 . Il988l) and they extend in mass only up to 12OM , although Georges Meynet 
(2010, private communication) says that there is no reason why higher mass models 
couldn't be made available if there was interest in these. What is the expected lifetime 
of a 300M Q st ar? As one approaches the Eddington limit, L ~ m rather than some 
steeper power ( Owocki & van Marlell2008h . and since the main-sequence lifetime r ms 
is going to go something like r ms ~ m/L we would expect the lifetimes to be fairly 
invariant with m ass. At 80-120M o t ~ m 5 . The lifetime of a 120M o star is 3.1 Myr at 
solar metallicity (Me vnet & Maededl2003l) . So, one would expect a lifetime somewhere 
between around 2 Myr for a 300M o star. 

Examination of Table 1 reveals that there are very few clusters and associations 
we know of that are that young. So, the fi rst thing to con s ider i s that if one is going to 
perform a statistical test along the lines of lOey & Clar ke ( 2005) then one really has to 
restrict one's self only to these few regions that are young enough. 
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The second point I wish to make is that we really don't know the masses of 
the highest mass s tars in these associations very well. Take the example of R136al. 
iMassev & Hunted (1 19981) makes a very conservative estimate of 150M G for its mass, 
but this required again extrapolating the mass-luminosity relationship further than the 
models go. As they note, the effective temperatures of such "Of stars on steroids" 
(which exhibit hydrogen-rich WN-type emission spectra) are very poorly known, and 
without a good effective temperature you don't know the bolometric luminosity very 
well, and without the bolometric luminosity you can't get the mass even by extrapolat- 
ing the mass-luminosity relationship. Since the time of the conference, ICrowther et al.l 
(2010) has made a new effective temperature determination by obtaining VLT optical 
data and fitting it (along with archival HST data). Th ey obtain a much hig her tem- 
perature (53,000 K) than that found in the analysis by Ide Koter et"ai1 (Il997l) (45.000 



K), al ong with a correspondingly higher luminosity and higher mass. ICrowther et al. 



( 20101) now estimate the (initial) mass of R136al be 300M G . I was tickled when I first 
saw the paper because during my talk I noted that one could argue for a higher tem- 
perature for this star and that indeed that a modern analysis with CMFGEN might well 
reveal it to be a 300M o and not 150M o . Bazinga! So, this is a very well done and 
interesting study, but when interviewed by the press I got to say that I wasn't surprised 
by it. 

Even so, this is har dly the final word on the subject. The stars in the core of R136 
are quite crowded, and dde Koter et alJH997l) argue that their optical HST spectrum is 
contaminated by neighboring stars. It is doubtful that the ground based spectrum, ob- 
tained even with AO, is without contamination. De Koter et al. (1997) also present 
strong arguments as to why the effective temperatures must be less than 47500 K. But 
I'll let the stellar atmosphere pundits slug that one out.) How significant this contami- 
nation may be, no one knows yet. Similarly, there is always the possibility that this star 
is an undetected binary. But, even so, it pretty well puts to rest the 150M o upper mass 
limit myth that has existed for the past few years. 

The final point to make about the upper mass limit is to remind the reader that age 
can do a very nice job truncating a luminosity function — the truncation does not h ave to 
be due to an upper mass cutoff. The Arches cluster provides an example of this. iFigerl 
(2005) cited the truncation of the K-band luminosity function to be proof that there was 
an upper mass cutoff of about 150M G : extrapolating the IMF suggested there should 
be 18 stars with masses > 130M o while none were fo und. But the picture was greatly 
clarified thanks to spectroscopy. Martins et al.1 ((2008) obtained spectra of stars, placed 
them on an HRD and concluded that "All stars are 2-4 Myr old". So, one expects that 
the luminosity function (and mass function) will be truncated. The region is just a bit 
too old to shed any light on the upper mass limit. 

I will say that my prejudice is that there certainly is some upper mass limit. Some- 
thing must limit the mass of a massive star. Our colleagues in the high mass star for- 
mation business here at this conference may tell us what. 



4. Summary 

This conference has certainly been held in a beautiful place! But I think Sedona also 
serves as a warning to us: having a number of pieces of indirect, weak evidence doesn't 
prove the variability of the upper end of the IMF or that 150M G is the upper mass limit, 
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any more than anecdotal stories prove the existence of UFOs and cosmic vortexes. They 
may be fun to believe in, but one has to look at the evidence critically. 

Acknowledgments. I'm grateful to Deidre Hunter and Michael Meyer for critical 
comments on the talk and manuscript, and for many useful discussions over the years. 
I'm also grateful to the conference organizers for bringing us all together here for such 
a pleasant meeting. 
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